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Near-Infrared Spectroscopy of Very Low-Luminosity Young Stellar 
Objects in the Taurus Molecular Cloud 
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Abstract 

We have carried out near-infrared spectroscopic observations of 23 very low-luminosity young stellar 
object (YSO) candidates and 5 their companions in Heiles Cloud 2, one of the densest parts of the Taurus 
molecular cloud. Twelve objects were confirmed as YSOs by Br7 feature. The effective temperatures of the 
YSOs and of the companions are estimated from the 2.26 //m feature, the 2.21 /jm feature, and the H2O 
band strengths. Detailed comparisons of our photometric and spectroscopic observations with evolutionary 
tracks on the HR diagram suggest some objects to be very low-mass YSOs. 
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1. Introduction 

Recent optical and near-infrared photometric 
studies have revealed faint populations of young 
stellar object (YSO) candidates in low-mass star 
forming regions (j^00 [cite]cite.ComeronComeron ct 
al.(1993); 000 [cite]cite.Strom95Strom et al.(1995); 
300 [cite]cite.ITGItoh et al.(1996), hereafter ITG; 
(]00 [cite] cite. BarsonyBarsony et al.(1997); ©00 
[cite]cite.Oasa990asa et al.(1999)), in intermediate- 
mass star forming regions (©00 [cite] cite. AspinAspin 
ct al.(1994)), and even in high-mass star form- 
ing regions (©00 [cite]cite.KaifuKaifu ct al.(2000); 
000 [cite]cite.Lucas00Lucas & Roche(2000); ©00 
[cite]cite.Oasa02Oasa et al(2002)). Their faintness 
is suggestive of low-mass; they may be very low-mass 
young stars near the stellar/substellar boundary, young 
brown dwarfs, or even free floating planets. 

From the photometric observations alone, how- 
ever, it is impossible to simultaneously determine 
the mass and age of a YSO. Near-infrared spec- 
troscopy of YSO candidates is necessary to over- 
come this difficulty (©00 [cite] cite. Greene95Greene 
& Meyer (1995); ©00 [cite]cite.Luhman97Luhman 
et al.(1997); ©00 [cite] cite. LuhmanRiekeLuhman 
& Rieke(1998); |(J)0 [cite] cite. Luhman98Luhman et 
al.(1998); ©00 [cite] cite. WilkingWilking et al.(1999); 
000 [cite]cite.CushingCushing ct al.(2000); ©00 
[citc]cite.Lucas01Lucas et al.(2001)). Spectra of some 
faint YSO candidates exhibit the absorption features 
distinctive to late spectral type, implying young brown 



dwarfs. While such very low-mass objects may be 
ubiquitous in star-forming regions, detail of the formation 
process of such objects, for example mass function, object 
den sity, and disk property, is still unknown. 

te]cite.Briceno98Briceno (1998) ([cite]cite.Briceno98Briceno 
(1998)) have carried out an optical search for very low- 
mass YSOs in the L1495, L1529, L1551, and B209 
regions in the Taurus molecular cloud. From photom- 
etry with spectroscopy, they found 9 new YSOs in 
the clouds, half of them have very late spectral types, 
implying very low-mass objects (0.05 M Q - 0.25 Mq). 
§e]cite.Luhman00Luhman(2000) ([cite]cite.Luhman00Luhman(2000)) 
further investigated very low-mass YSOs in the same 
clouds by combining optical imaging with near-infrared 
spectroscopy. They found that the mass function of 
the YSOs in these regions has a peak at 0.8 Af Q and 
is relatively flat between 0.1 M and 0.8 Af Q range. 
In these papers, however, because the targets for the 
spectroscopy were fully or partially selected based on the 
optical color-magnitude diagram, the sample may not be 
complete especially for embedded very low-mass objects. 

ITG conducted a near-infrared survey of the central 1° x 
1° region of Heiles Cloud 2 in the Taurus molecular cloud, 
one of the best-studied low-mass star forming regions, 
with a limiting magnitude of 13.4 mag in the if -band. 
Fifty YSO candidates were identified by their intrinsic 
red color on the (J — H , H — K) color-color diagram, 
following the scheme discussed by |t]e]cite.Strom93Strom 
ct al.(1993) ([cite]cite.Strom93Strom ct al.(1993)). 
Successive high-resolution imaging survey discovered 
5 companion candidates around the YSO candidates 



2 



Itoh et al. 



(§00 [cite]cite.ITNItoh et al.(1999)). Faintness of some 
YSOs and their companion may be an indicative of the 
low-mass of the objects. 

We describe here near-infrared spectroscopic follow-up 
of these faint YSO candidates in the Heiles Cloud 2. The 
observations are described in §2, and data reduction pro- 
cedures in §3. In §4, we derive effective temperatures of 
the YSOs mainly from the 2.21 //m feature, the 2.26 /tni 
feature, and the H2O absorption band. We also calculate 
photospheric luminosity of the YSOs from the previous 
photometry, then plot the YSOs on the HR diagram. 

2. Observations 

2.1. UKIRT observations 

Twenty-one YSO candidates were selected for the spec- 
troscopic observations (table |IJ). The sample is nearly 
flux-limited; The A"-band magnitudes range from 8 to 12, 
up to 4 mag fainter than that of classical T Tauri Stars 
(CTTSs) in the Taurus molecular cloud. We missed ITG 
8 and ITG 34 probably due to errors in the coordinates. 
As template stars, 3 late-type CTTSs, 1 class I source, 
and 7 late-type dwarfs were also observed. 

Spectroscopic observations have been carried out dur- 
ing 1996 November 24 - 26 using the UK Infrared 
Telescope (UKIRT) at the summit of Mauna Kea, 
with the Cooled Grating Spectrometer 4 (CGS4, ©00 
[cite]cite.MountainMountain et al.(1991)). CGS4 has a 
256x256 InSb array with a spatial scale of l."22 pixel" 1 . 
The 75 line mm" 1 grating was used with the l."2 slit, 
providing wavelength coverage of 1.8 /im - 2.5 /im with 
resolution R(~ A/AA) of 300 at 2.21 /im. The integration 
time used for each source was typically 1-20 seconds de- 
pending on source brightness. In most cases, 96 exposures 
were taken for each object. Nodding of the telescope was 
carried out approximately 16 arcsec along the slit for sky 
subtraction. For the YSO candidates and CTTSs, SAO 
76542 (A2V) was observed for correction of the effects 
of telluric absorption. For late-type dwarfs, we observed 
stars with spectral types of A0V - A3V at similar air- 
mass. Exposures of an incandescent lamp on and off were 
taken at the start of each night as dome flats. Exposures 
of a krypton lamp were taken every three or four hours 
for wavelength calibrations. 

2.2. Subaru observations 

The observed objects are 5 binary candidates listed 
in §e]cite.ITNItoh et al.(1999) ([cite]cite.ITNItoh et 
al.(1999)) (table 0). The A-band magnitudes of the com- 
panions range from 13 to 16. Six late- type dwarfs were 
also observed as templates, among which 3 latest dwarfs 
were observed both with UKIRT and Subaru. 

Spectroscopic observations have been carried out on 
2000 December 4-5 with the Infrared Camera and 
Spectrograph (IRCS, §00 [cite]cite.Tokunaga98Tokunaga 
et al.(1998); §30 [cite]cite.KobayashiKobayashi et 
al.(2000)) on the Subaru Telescope at the summit of 
Mauna Kea, Hawaii. IRCS has a 1024x1024 InSb array 
with spatial scale of 0."058 pixel -1 . Typical seeing size 



was 0"4 with a stable condition for both observing dates, 
so that all binaries were well separated. The grating 
provided a wavelength coverage of 2.0 /im - 2.5 /im. The 
resolving power R was around 350 at 2.2 /tm. The slit 
width we used is 0"6, slightly wider than the seeing size, 
so that the effective spectral resolution might change 
by seeing size. We measured FWHMs of the absorption 
features in each spectrum for each object, and confirmed 
that the FWHMs do not change in 1 pixel resolution 
(~ 6 A). Therefore, the effective spectral resolution was 
stable during observations of each target. The integration 
time used for each source was typically 60 - 300 seconds 
depending on source brightness. 8 or 12 exposures 
were taken for each object with the telescope dithered 
approximately 7"5 along the slit for sky subtraction. 
For the binary candidates, SAO 76542 was observed 
for correction of the effects of telluric absorption. For 
late-type dwarfs, we observed stars with spectral types 
of A0V - A3V at similar airmass. Exposures of an 
incandescent lamp on and off were taken at the end of 
each night as dome flats. Exposures of an argon lamp 
were taken for wavelength calibrations at the end of each 
night. 

3. Data Reduction and Results 

The Image Reduction and Analysis Facility (IRAF) 
software was used for all data reduction. Reduction pro- 
cedure is similar for both UKIRT and Subaru data. First, 
a dithered pair of object frames were subtracted by each 
other, then divided by flat fields. Next, the data frames 
were geometrically transformed to correct the curvature 
of the slit image caused by the grating. The solutions of 
the geometric transformation and the wavelength calibra- 
tion were derived from the spectrum of a krypton lamp or 
argon lamp taken closest in time to the object. Then, in- 
dividual spectra were extracted from the transformed im- 
ages using the APALL task. The region where the inten- 
sity of the object was more than 20 % of the peak intensity 
at each wavelength was summed into a one-dimensional 
spectrum. Extracted spectra were then normalized and 
combined to produce the final spectra. Low signal-to- 
noise spectra due to the tracking error of the telescope 
were rejected. The object spectrum was divided by the 
standard star spectrum, in which Brackett7 absorption 
line at 2.166 /tm was removed by interpolating across the 
adjacent continuum with the SPLOT task. Finally the 
spectrum was multiplied by a blackbody spectrum of the 
temperature appropriate to the spectral type of the stan- 
dard star (§00 [cite]cite.TokunagaTokunaga(2000)). 

The A-band spectra of the 23 YSO candidates, 5 com- 
panion candidates, 3 CTTSs, a protostar, and 10 late-type 
dwarfs are shown in figures [l] and §. Prominent features 
in the spectra are the HI Br7 (2.17 /im), 2.21 /tm feature, 
2.26 /tm feature, CO band (2.29 /tm and longer), and H2O 
absorption (<2.15 /tm and >2.3 /tm). 

Equivalent widths of the Br7 line, the 2.21 /tm feature, 
and the 2.26 /tm feature were measured with SPLOT task 
by Gaussian fitting. On the other hand, equivalent widths 
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Table 1. The YSO sample in Heiles Cloud 2 
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ITG No. g(1950) (5(1950) K J-H H-K Av Obs.* Identification^ 



2 


4 h 


34 m 57 s 


25°53'01" 


10 


,05±0 


.01 





,84d=0. 


01 





60±0, 


01 


0.14±0. 


49 


U 






4 


4 


35 


15.0 


2G 


01 


28 


11 


,05±0 


.02 


2 


07±0 


.05 


1 


,29zbO, 


,03 


11.42±0 


.55 


U 






5 


4 


35 


17.1 


25 


4G 


06 


8. 


31±0 


.00 


1 


,17±0 


.01 





.72±0. 


.01 


3.38±0. 


,47 


u 


Kim 1-19 




6 


4 


35 


17.5 


26 


03 


19 


10 


,37±0 


,01 





96±0 


.02 





,87±0 


.02 


0.14±0 


,49 


u 


GM Tau 




9A 


4 


35 


56.6 


25 


27 


17 


12 


,53±0 


.03 





.88±0. 


,04 





66±0 


,04 


0.27zb0. 


,42 


s 






9B 


4 


35 


56.6 


25 


27 


21 


14 


,45±0. 


.03 





70±0 


.04 





,70±0 


,04 


0.00±0 


.39 


s 


separation = 


4."31 


13 


4 


36 


15.2 


25 


31 


44 


11 


69±0 


.02 





.84±0. 


,03 





,55±0 


,03 


0.36±0 


.51 


u 






15A 


4 


3G 


40.5 


25 


56 


02 


9 


02±0 


.00 





.86±0. 


,01 





,62dz0 


,01 


0.19±0. 


.47 


u,s 






15B 


4 


36 


40.4 


25 


56 


06 


14 


50±0 


.01 


1 


OldzO, 


,02 





,73dz0 


,02 


1.48±0. 


.39 


s 


separation = 


2."99 


17 


4 


3G 


43.1 


25 


55 


49 


10 


,37±0 


.01 





96d=0 


,02 





.82±0. 


,02 


0.45zb0. 


.49 


u 






18 


4 


36 


46.3 


25 


41 


02 


10 


,27±0 


.01 


1 


96±0 


,03 


1 


,18±0, 


,02 


10.6±0 


.50 


u 






21 


4 


36 


57.5 


25 


50 


38 


10 


,83±0 


.01 


1 


,49dzO 


,03 





94±0 


,02 


6.21±0 


.51 


u 


GKH 5 




24 


4 


37 


03.6 


26 


01 


2G 


13 


,05±0 


.05 


1 


,01±0. 


,09 





,70±0. 


,09 


1.65±0 


.82 


u 






25A 


4 


37 


03.8 


25 


59 


38 


8 


,97±0 


.00 


1 


97±0 


,01 


1 


37±0 


,01 


9.82dz0 


.49 


u,s 


IRAS 043704 


-2559 


25B 


4 


37 


03.6 


25 


59 


36 


13 


,38±0 


.01 


2 


14±0 


,01 


1 


.28±0. 


,01 


12.29d=0 


.38 


s 


separation = 


4f 29 


27 


4 


37 


32.7 


25 


52 


34 


8 


,46±0 


.00 


1 


,04±0. 


,00 





64dz0 


,00 


2.28±0 


.47 


u 






28 


4 


37 


46.2 


25 


13 


16 


11 


,87±0 


.02 





,81±0. 


,04 





53±0 


,04 


0.09dz0 


.54 


u 






29 


4 


3 1 


52.5 


25 


16 


19 


10 


,74±0 


.01 





.86±0. 


,01 





,54±0. 


,01 


0.63dz0 


.49 


u 






33A 


4 


38 


04.2 


25 


50 


23 


11 


,50±0 


.02 


1 


A 1 I A 

,41±0 


04 


1 


19±0 


no 

.03 


3.94±0 


,55 


TT C 

U,b 






33B 


4 


38 


04.1 


25 


50 


28 


16 


,05±0 


.02 





79±0 


,04 





,54d=0. 


.03 


0.00±0 


,41 


s 


separation = 


5."17 


36 


4 


38 


10.0 


25 


12 


27 


11 


.OldbO 


.02 





86±0. 


,02 





,54±0 


,02 


0.63±0 


,49 


u 






39 


4 


38 


16.4 


25 


28 


24 


11 


,80±0 


.02 


1 


,01d=0. 


,03 





66zb0 


,03 


1.88±0. 


,51 


u 






40 


4 


38 


20.8 


25 


38 


10 


11 


,47±0 


.02 


2 


.80±0. 


,18 


1 


,99zb0 


,05 


16.68±0 


,83 


u 


GKH 32 




41 


4 


38 


34.6 


25 


50 


4G 


9 


,72±0 


.01 


1 


66dz0 


,01 


1 


llzbO, 


,01 


7.38±0 


,49 


u 


IRAS 043854 


-2550 


43 


4 


38 


39.8 


25 


27 


30 


11 


,24±0 


.01 





.84±0. 


,01 





53±0 


,01 


2.97zb0. 


,55 


u 






45A 


4 


38 


45.3 


25 


42 


37 


12 


,77±0 


.05 





90±0 


,06 





73±0 


,07 


0.21zb0. 


,52 


s 






45B 


4 


38 


45.2 


25 


42 


37 


15 


,76±0 


.05 


1 


33±0 


,08 





,91±0 


,09 


4.42±0. 


,62 


s 


separation = 


2'.'29 


46 


4 


38 


52.3 


25 


47 


16 


11 


.OldbO 


.01 





,83±0 


,01 





,52±0 


,01 


0.41±0. 


49 


u 







* U: UKIRT, S: Subaru 

t GKH: Bs]cite.GomezGomez et al.(1994) ([cite] cite. GomezGomez et al.(1994)), Kim: |e]cite.KimKim(1990) ([cite]cite.KimKim(1990)) 



equivalent widths and band strengths are within observa- 
tional uncertainties. 

4. Discussion 

4-1- Identification as YSOs using the Brj feature 

We identified 12 sources as YSOs, whose Br7 feature is 
in emission or in flat (featureless). We regard a feature- 
less spectrum around the Br7 line with uncertainty less 
than 2 A in equivalent width as flat. The line emission is 
probably due to accretion of matter from the circumstel- 
lar disk onto the star (e.g. [)]00 [citejcite.NajitaNajita et 
al.(1996)). 

The other 11 objects with Br7 absorption are, on 
the other hand, likely early- type field-stars. Equivalent 
widths of the Br7 absorption feature are 4 A- 15 A, 
consistent with the spectral types of B, A, and early F 
(§00 [cite] cite. AliAli et al.(1995)). Their locations on 
the color-color diagram are close to that of early-type 
field stars. Figure § shows the color-color diagram of 
the YSO candidates and field stars identified by ITG 
toward Heiles Cloud 2. The objects with the Br7 ab- 
sorption tend to be located near the boundary between 



of the CO(2-0) band and the CO(4-2) band were calcu- 
lated by simple integration of the absorption intensity. 
The uncertainties were estimated from the continuum fit, 
in which locating the continuum level was the main factor 
contributing to the equivalent width uncertainty. We 
also calculated reddening- independent indices of the H2O 
band strength Q, following |t]e]cite.WilkingWilking 
et al.(1999) ([cite] cite. WilkingWilking et 

al.(1999)). For Koornneef's extinction law (§00 
[cite]citc.KoornnecfKoornncef(1983)), the Q index is 
written as 

where F\, F2, and F3 are flux densities between 2.07 /xm 
- 2.13 fim, 2.267 fim - 2.285 fim, and 2.40 ^m - 2.45 /zm, 
respectively. 

The measured equivalent widths of the features and 
the strengths of the bands are tabulated in table |2|. The 
signal-to-noise ratio derived from deviations between each 
exposure is typically 60 (table §). For the latest dwarfs 
taken both with UKIRT and Subaru, shapes of the spectra 
are similar for both observations and most of the measured 
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Table 2. Equivalent widths of the features and strengths of the band of the YSO candidates, the CTTSs, and the M dwarfs. 



Object 


Br 7 


2.21 /.an 


2.26 fim 


CO(2-0) 


CO(4-2) 


Q 


Obs.* 


S/N 




[A] 


[A] 


[A] 

L J 


[A] 

L J 


[A] 

L J 








ITG 2 


0.0±0.8 


2.18±0.39 


1.00+0.16 


6.20 


7.80 


0.51+0.00 


u 


70 


ITG 4 


6.80±0.25 


0.0±0.3 


0.0+0.1 






0.98+0.01 


u 


90 


ITG 5 


4.24±1.50 


0.80±0.26 


1.44+0.28 


2.71 


1.97 


0.83+0.01 


u 


90 


ITG 6 


-2.25±0.62 


0.98±0.23 


0.53+0.22 


2.30 


3.79 


0.81+0.01 


u 


90 


ITG 9A 


0.0±0.5 


2.42±0.04 


2.38+0.08 


3.45 


4.72 


0.89+0.02 


s 


300 


ITG 9B 


0±3 


0±3 


0.0+0.2 


2.71 


2.35 


0.88+0.09 


s 


60 


ITG 13 


14.08±0.70 


0.0±0.3 


0.0+0.4 






0.94+0.01 


u 


60 


ITG 15A&B 


-0.79±0.79 


3.29±1.00 


4.12+1.24 


6.73 


10.39 


0.64+0.01 


u 


80 


ITG 15A 


-3.93±1.08 


3.29±0.24 


3.19+0.30 


5.88 


9.79 


0.61+0.00 


s 


460 


ITG 15B 


0.0±2.5 


3 53+2 9? 

'.J • ' ) * } _i_ _.-/_■ 


14 19+12 37 


25.90 


26.40 


1 65+nti 

1 - UJ -0.82 


s 


9 


ITG 17 


0.0±0.5 


1 96±0 37 


2.48+0.17 


6.79 


10.02 


57+0 01 


TJ 


30 


ITG 18 


8.25±0.62 


0.0±0.3 


0.0+0.3 






98+0 01 


TJ 


80 


ITG 21 


-3.8±3.8 


2 22+0 78 


2 31+0 44 

_ . ' 7 J. 1 \ 1 ■ 11 


5.11 


8.50 


54+0 01 

\J . t7^t_LU . \J ±. 


TJ 


70 


ITG 24 


6 60±1 69 


0.0+0.5 


58+0 1 2 






1 03+0 01 


TJ 


40 


ITG 25A&B 


-5 34±0 91 


1 94+0 51 


2 54+0 50 


3.72 


5.43 


87+0 03 


TJ 


100 


ITG 25A 


-12.63±0.22 


3 38+0 16 


2.26+0.08 


3.52 


8.26 


89+0 02 


s 


290 


ITG 25B 


0.0±3.4 


1 41+0 60 


2 53+1 13 


9.61 


21.31 


73+0 09 


s 


60 


ITG 27 


0.0±1.2 


1 36+0 64 


1 33+0 70 


7.87 


8.50 


86+0 01 


TJ 


100 


ITG 28 


5.52±0.72 


0.75+0.15 


0.0+0.4 






1.00+0.01 


TJ 


50 


ITG 29 


11.54±0.80 


0.47+0.12 


0.0+0.1 






0.96+0.01 


TJ 


65 


ITG 33A&B 


-2.83±1.05 


3.68+0.43 


3.76+0.37 


6.16 


8.56 


0.74+0.01 


U 


55 


ITG 33A 


-5.85±0.85 


3.78+0.54 


3.92+0.47 


6.02 


7.35 


0.73+0.01 


s 


660 


ITG 33B 


0.0±2.9 


7.54+4.90 




0+9 




0.90+0.04 


s 


20 


ITG 36 


5.14±1.63 


1.20+0.43 


1.26+0.33 






0.91+0.01 


u 


50 


ITG 39 


7.40±0.70 


0.0+0.2 


0.0+0.2 






1.00+0.02 


u 


60 


ITG 40 


0.0±0.9 


4.15+1.00 


1.78+0.53 


4.95 


10.76 


0.67+0.02 


u 


30 


ITG 41 


0.0+1.0 


2.31+0.50 


2.08+0.56 


3.43 


6.53 


0.93+0.01 


u 


100 


ITG 43 


7.68±0.98 


0.0+0.4 


0.0+0.2 






0.97+0.01 


u 


60 


ITG 45A 


8.26±0.96 


0.0+0.88 


0.0+1.8 






0.88+0.07 


s 


80 


ITG 45B 


0±25 


0+4 


0+9 






^^ — 0.83 


s 


6 


ITG 46 


0.0+1.0 


1.46+0.50 


1.34+0.33 


11.13 


13.69 


0.80+0.01 


TJ 


70 


DD Tau (Ml) 


-2.53±0.35 


1.99+0.42 


1.84+0.25 


3.21 


4.93 


0.82+0.01 


TJ 


100 


GH Tau (M2) 


-0.58±0.57 


2.91+0.46 


2.81+0.48 


5.01 
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TJ 
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-1.02±0.46 
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TJ 
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0.39+0.07 


14.32 


13.54 
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Fig. 1. (a) ii"-band spectra of the YSO candidates taken with Fig. 1. (c) i\"-band spectra of the YSO candidates taken with 

UKIRT. The spectra are normalized by the flux between 2.18 UKIRT. The spectra are offset in steps of 0.5. 

(im and 2.20 fim, then offset in steps of 0.5. 




Wavelength (f^m) 



Fig. 1. (b) X-band spectra of the YSO candidates taken 
with UKIRT. The spectra are offset in steps of 0.5. 



Fig. 1. (d) i\"-band spectra of the YSO binary candidates 
taken with Subaru. Additive constants for the spectra are 0, 
1, 1.5, and 2, respectively. 
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Wavelength (j-im) Wavelength (/i.m) 

Fig. 1. (e) _ff-band spectra of the YSO binary candidates Fig. 2. (b) The i^-band spectra of late-type dwarfs taken 

taken with Subaru. Additive constants are 0, 1, 1.5, 2, 2.5, with UKIRT. The spectra are offset in steps of 0.5. 
and 3. 



3 - 



Br-/ 3.21 2.26 CO - 
H-0 H„0 



I RAS04 365 + 2535 




2. i 



2.2 2.3 
Wavelength (/im) 



2. -I 



Fig. 2. (a) A*-band spectra of CTTSs and a protostar taken 
with UKIRT. The spectra are normalized by the flux between 
2.18 fim and 2.20 /Ltm. Additive constants for the spectra are 
0, 0.5, 1, and 2, respectively. 
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Fig. 2. (c) i^-band spectra of late-type dwarfs taken with 
Subaru. The spectra are offset in steps of 0.5. 
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□ Bry absorption 

x No spectrum available 
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the "YSO region" and the "field-star region". Therefore, 
these objects may be early-type field-stars and have been 
misidentified as YSO candidates by photometric classifi- 
cation. The star count model of |e] cite. Jones81 Jones et 
al.(1981) ([cite]cite.Jones81 Jones et al.(1981)) predicts 9 
B- and A-type stars for if < 12 in a 1° x 1° region to- 
ward Heiles Cloud 2. This predicted number is consistent 
with the observational number of the early type stars. 
Note that the number of such early-type stars does not 
increase for K > 12 toward high latitude direction such as 
the Taurus molecular cloud. 

The uncertainty of the slope of the reddening vector 
may also account for the misidentification. For the Taurus 
molecular cloud, the smallest value of E(J — H)/E(H — K) 
is 1.58 (§00 [cite]cite.EliasElias(1978)), whereas the high- 
est is 2.0 (§00 [cite] cite. Gomez Gomez et al.(1994)). From 
photometric study alone, identification of YSOs depends 
on the extinction law especially for the objects near the 
boundary between the YSO region and the field star re- 
gion on the color-color diagram. 

4-2. Effective Temperature 

We first describe the possible if-band features with 
which the effective temperatures of the M dwarfs and gi- 
ants are derived. We then apply this method to the YSOs 
to derive their effective temperatures. 

For cooler objects the conversion from MK spectral 
type to effective temperature is not straightforward 
because models do not exist to give the spectral 
type vs. effective temperature for pre-main sequence 
stars. Because the YSOs have luminosity class be- 
tween dwarfs and giants, as we discuss below, we use 
the conversions both for dwarfs and for giants. In 
this paper we have adopted the conversions given by 

|e]cite.TokunagaTokunaga(2000) ([citeJcite.TokunagaTokunaga^QflSOj^ccts, and the CTTSs whose ncar-infrarcd excess is 
both for dwarfs and giants, small. In the "D" (disk) region, a part of CTTSs are located 

|elcite.Bessell91Bessell(1991) ([citelcite.Bessell91Bessell(1991)) Protostars are located in the " E " ( envel °P e ) re g ion ° n thc 
for late M dwarfs, and |tp]cite.FluksFluks et 
al.(1994) ([citc]citc.FluksFluks et al.(1994)) for late 
M giants. As can be seen in figure |[ the difference in 
effective temperature is as much as 800K between M 
dwarfs and M giants even for the same spectral type. 



Fig. 3. The color-color diagram for Hcilcs Cloud 2. 
The intrinsic colors of main-sequence stars, giants ( ] 30 
[cite] cite. Bessell88Bessell & Brett(1988)), supergiants ('30 
[cite]cite.TokunagaTokunaga(2000)), and CTTSs ('30 
[cite]cite.Meyer97Meyer et al.(1997)) are plotted with thc 
reddening vector (MOO [cite]cite.KoornneefKoornnccf(1983); 
E(J-H)/E(H-K)=1.7 in the CIT system). All colors 
are transformed to the Johnson/Glass system (HOO 
[cite] cite. Bessell88Bessell & Brett(1988)). Following 
the scheme discussed by me]cite.Strom93Strom et 
al.(1993) ([cite] cite. Strom93Strom et al.(1993)), thc 
near-infrared sources are classified into three groups. In 
the "P" (photosphere) region on the color-color diagram, 
located are main sequence stars, giants, weak-line T Tauri 



4.2.1. 2.21 /Ltm feature 

The 2.21 /im feature consists mainly of Na, Sc, 
Si, and V absorption lines for late-type objects (§00 
[cite] cite. RamirezRamircz et al.(1997)). The low exci- 
tation energy of these lines makes it a strong absorp- 
tion feature in M type stars. Figure [5] shows equi- 
valent widths of the 2.21 /im feature of dwarfs and gi- 
ants as a function of effective temperature. The equi- 
valent widths increase with decreasing effective tem- 
perature. Moreover, dwarfs have larger equivalent 
widths than giants at the same effective tempera- 
ture. The 2.21 fim feature has, therefore, a strong 
dependence on effective temperature and weak depen- 
dence on gravity (§00 [cite]cite.KleinmannKleinmann & 
Hall(1986); §00 [citc]citc.TerndrupTerndrup et al.(1991); 
§00 [cite] cite. RamirezRamirez et al.(1997)). 



color-color diagram. The objects located in the "D" region 
or the "E" region were identified as YSO candidates by the 
previous photometric study (ITG). The objects without 
absorption (i.e. with emission or in flat) in the Br7 line 
are denoted by filled squares, as well as the objects with 
absorption by open squares. The YSO candidates without 
the spectroscopic observations are plotted by crosses. Field 
stars identified by ITG arc shown by dots. 
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Fig. 4. Relation between effective temperature and 
spectral type for dwarfs and giants. Data are taken from 
I s]cite.TokunagaTokunaga(2000) ([cite]cite.TokunagaTokunaga(2 
|c]citc.WilkingWilking et al.(1999) ([cite]cite.WilkingWilking 
et al.(1999)), HcJcite.LuhmanRiekeLuhman & 

Rieke(1998) ([cite]cite.LuhmanRiekeLuhman & Rieke(1998)), 
[c]citc.Bessell91Bessell(1991) ([cite]cite.Bessell91Bessell(1991)), 
fsjcitc.DyckDyck et al.(1998) ([cite]cite.DyckDyck 
ct al.(1998)), ^cite.PerrinPerrin et 

^,1.(1998) ([cite]cite.PerrinPerrin et al.(1998)), and 
|e]citc.FluksFluks et al.(1994) ([cite]cite.FluksFluks ct 
al.(1994)). 

Near-infrared Na lines of CTTSs have symmet- 
ric profile, indicating photospheric origin (^00 
[cite]cite.Greene97Greene & Lada(f997)), while op- 
tical Na lines are red-shifted, indicating outflow origin 
(@00 [cite]cite.GuUbringGuUbring et al.(1996)). 

4.2.2. 2.26 fj,m feature 

The 2.26 /mi feature consists mainly of Ca, Si, and Ti 
lines for late-type objects (^00 [citeJcite.RamirezRamirez 
ct al.(1997)). Figure ^| shows equivalent widths of the 
2.26 /im feature of dwarfs and giants as a function of ef- 
fective temperature. The 2.26 /im feature has a peak in its 
equivalent width around 3500 K, then its width decreases 
with decreasing temperature, due to moderate excitation 
energy of Ca. This feature has also strong dependence on 
effective temperature, and weak dependence on gravity 
(^|00 [cite] cite. Ramirez Ramirez et al.(1997)). 

4-2.3. Ratio of Metallic features 

Use of ratio of feature strengths to derive effective 
temperature can avoid not only the effect of metal- 
licity but also the effect of veiling if features whose 
wavelengths are close with each other are selected (^00 
[cite]cite.Meyer96Meyer(1996)). We have checked if the 
2.21 fim feature and the 2.26 /im feature can be used as 



Fig. 5. Equivalent widths of the 2.21 fim fea- 
ture as a function of effective temperature. Data 
0)),are taken from |I]e]cite.LuhmanRiekeLuhman & 
Rickc(1998) ([citcJcitc.LuhmanRiekeLuhman & 

Rieke(1998)), |e]cite.AliAli ct al.(1995) ([citc]citc.AliAli 
et al.(1995)), HcJcite.KleinmannKlcinmann 

& Hall(1986) X[citc]citc.KlcinmannKlcinmarm 

& Hall(1986)), rp]citc.RamirezRamirez ct 

al.(1997) ([citeJcite.RamirezRamirez et al.(1997)), and 
this paper. Filled symbols represent the equivalent widths 
for dwarfs and open symbols for giants. A solid line 
shows a fitted line for dwarfs with T e g > 2300 K, and a 
dotted line for giants listed in ysJcite.RamirezRamircz ct 
al.(1997) ([citeJcite.RamirezRamirez ct al.(1997)). 

a ratio index of effective temperature. Figure Q shows 
the 2.26 fim / 2.21 /im feature ratio for dwarfs and gi- 
ants as a function of effective temperature. Although the 
scattering is relatively large and no data is available for 
cool giants, the ratio decreases as the temperature de- 
creases. Dwarfs and giants have similar values, so that 
surface gravity does not appear to significantly affect this 
feature ratio. Therefore, this feature ratio is a sensitive 
indicator of effective temperature for T e ff < 3300 K. For 
2.26 /im / 2.21 /tm > 1, only lower limits of effective tem- 
perature can be estimated. 

4.2.4. H 2 

The H2O absorption bands in the near-infrared wave- 
lengths are a sensitive function of effective temper- 
ature for cool stars ([]00 [cite]cite.AaronsonAaronson 
et al.(1978); ©00 [citc]citc.Jones94Jones et al.(1994)). 
Figure || shows the water index Q of the dwarfs 
taken from our observations, as well as those of 
dwarfs and giants taken from [t]e]cite.LanconLangon 
& Rocca-Volmerange(1992) ([cite] cite. LanconLangon & 
Rocca-Volmerange(1992)). A dotted line represents 
a fitted line of the Q index for the dwarfs and 
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Fig. 6. Equivalent width of the 2.26 fim feature as a func- 
tion of effective temperature. Filled symbols represent 
the equivalent widths for dwarfs and open symbols for 
giants. Solid lines show fitted lines for dwarfs, and a 
dotted line for giants listed in |I|c]cite.RamirezRamirez et 
al.(1997) ([cite]cite.RamirezRamirez et al.(1997)). The inten- 
sity has a peak at 3500 K for dwarfs. 

a dashed line for the giants. In this figure, the 
relation between the Q index and spectral type 
for M dwarfs derived by |e]cite.WilkingWilking et 
al.(1999) ([cite]cite.WilkingWilking et al.(1999)) also 
shown by a solid line. While Wilking's relation is consis- 
tent with the plotted Q index for the dwarfs, the Q indices 
for giants are much smaller than the Wilking's relation. 
Since the YSOs have surface gravity between dwarfs and 
giants, as we discuss below, we should use both the re- 
lation for the dwarfs and that for the giants to estimate 
effective temperature of the YSOs. Note that the veiling 
effect by a circumstellar disk could increase the Q index, 
leading to higher temperatures. 

4-2.5. Effective Temperatures of the YSOs 

The effective temperatures of the observed YSOs and 
their companions are derived using the 2.26 /im / 2.21 
/an feature ratio and the Q index (T c ff(NIR) in table ||). 
The effective temperatures of most YSOs are less than 
4000 K, lower than that of typical CTTSs. 

The effective temperatures of the observed CTTSs are 
also derived from the near-infrared features. The effec- 
tive temperatures estimated from the near-infrared spec- 
troscopy are consistent with those converted from optical 
spectral types for all the CTTSs. 

4-3. Luminosity Class 

Since the 2.21 fj,m feature and the 2.26 /im fea- 
ture are primarily sensitive to effective temperature, and 



Fig. 7. Ratio of the 2.26 ^m feature to the 2.21 
fim feature as a function of effective tempera- 
ture. Data are taken from our UKIRT observations, 
Subaru observations, HcJcitc.LuhmanRickcLuhman 
&: Riekc(1998) ([citeJcite.LuhmanRiekeLuhman & 
Rieke(1998)), fle]cite.AKAli ct al.(1995) ([cite]cite.AliAli 
et al.(1995)), |I|c]cite.KleinmannKleinmann & 

Hall(1986) ([citcJcite.KleinmannKlcinmann &: 

Hall(1986)), and |t]e]cite.RamirezRamirez et 

al.(1997) ([cite]cite.RamirezRamirez et al.(1997)). Solid lines 
show the ratio for dwarfs, and a dotted line for giants, both 
derived from the fitted lines in figures |jj and ^. Dashed lines 
show the conversions used for estimate of the upper limit 
and the lower limit of effective temperature of the YSOs. 

the depth of the CO band varies with surface gravity, 
the luminosity class of the object is estimated from a 
plot of the equivalent widths of these atomic features 
against the widths of the CO band. Figure || shows 
the sum of the equivalent widths of the 2.21 /jm fea- 
ture and the 2.26 /im feature against the CO band equi- 
valent width of the YSOs, the CTTSs, and the late- 
type dwarfs. A sum of the equivalent widths of the CO 
(2-0) and CO (4-2) bands is used as a CO equivalent 
width. CO (3-1) band is not included because it con- 
tains a sodium line (|(]00 [cite]cite.KleinmannKleinmann 
& Hall(1986)). Relationships for giants and dwarfs (j^00 
[cite]cite.Greene95Greene & Meyer(1995)) are also shown 
in the figure by the solid lines. Their fits arc only applied 
to stars with spectral types < M6V and < M2III, and 
cannot be applied to very late-type stars. Points of very 
late-type dwarfs are close to the line for giants, since their 
equivalent width of the 2.26 /jm feature diminishes (see 
figure |). 

Most of the YSOs lie, like the CTTSs, between these 
two lines. Therefore most of the YSOs and the CTTSs 
have a luminosity class between dwarfs and giants, classi- 
fied as luminosity class IV. 
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Table 3. Effective temperatures and spectral types of the YSOs and the CTTSs. 



Object 
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T cff (Opt) # 


Sd (V)** 


Sd QUI 1 " 1 ' 


log Lbol" 






IAJ 


[Kl 


[Kl 




[K] 
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* T c ff estimated from the 2.26 jim / 2.21 (im feature ratio. 

t T e ff estimated from the the equivalent widths of the 2.21 /im feature. Because YSOs are subject to continuum emission and the strength 
of the line increases with lower temperature, only the upper limit of the effective temperature can be estimated. 

T e g estimated from the equivalent widths of the 2.26 /im feature. 
§ T c ff estimated from the Q index. The range of temperatures for T e g(Q) includes estimates for both dwarf and giant surface gravities. 
These effective temperatures may be upper limits, due to veiling effect. 
" T ef f estimated from the 2.26 /im / 2.21 (im feature ratio and the Q index. 

* T e ff converted from optical spectral type. 
** Estimated spectral type in the dwarf scale. 

tt Estimated spectral type in the giant scale. 

M Bolomctric luminosity for the midpoint of T c ff (NIR), estimated from the previous photometric study (see §1.4.1). 
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Fig. 8. The Q index as a function of effec- 
tive temperature. Data are taken from our 
observations and from jtjcJcite.LanconLangon & 
Rocca-Volmerange(1992) ([citc]cite.LanconLancpn & 
Rocca-Volmerange(1992)). A solid line shows the relation- 
ship between the Q indices and spectral type for dwarfs (HOO 
[citc]citc.WilkingWilking et al.(1999)). A dotted line shows 
a fitted line of the Q index for the dwarfs, and a dashed line 
for the giants. 
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Fig. 9. Luminosity class of the YSOs. The combined equi- 
valent width of the 2.21 fim feature and the 2.26 fim feature 
is plotted against the combined CO band equivalent width. 
Loci of dwarfs and giants ([jt)0 [cite]cite.Greene95Greene & 
Meyer(1995)) are shown by solid lines. 



4-4- Mass and Age of the YSOs 

Based on the effective temperature described above and 
the bolometric luminosity described below, the masses of 
the YSOs and their companions are estimated with recent 
evolutionary tracks on the HR diagram. 

4- 4-1. Bolometric Luminosity 

Bolometric luminosity of the YSOs is estimated from 
the J-band luminosity, because the J-band emission 
arises primarily from the photosphere of YSOs (j^OO 
[citc]citc.BcrtoutBcrtout ct al.(1988)). First, the J-band 
luminosity of the YSOs is corrected for interstellar extinc- 
tion on the color-color diagram (ITG). Next, three cases 
are assumed, where the fraction of the photospheric lumi- 
nosity to the total luminosity is 1.0, 0.8, and 0.5 in the 
J-band. Then, by comparing the photospheric luminosity 
of the YSOs with that of late-type dwarfs in the J-band 
(^00 [cite] cite. LeggettLeggctt ct al.(1996)), the bolometric 
hmiinosity of the YSO's photosphere is estimated 

This procedure has the following uncertainties. 



Extinction of each YSO 
on the color-color diagram. 



is 



corrected 
If the ex- 
tinction law of []e]cite.Bcsscll88Bessell 
& Brett(1988) ([cite]cite.Bessell88Bessell 
& Brctt(1988)) is used instead of 
|e] cite. KoornneefKoornneef (1983) ([cite]cite.KoornnecfKoornne 
a 6 % larger value is derived as visual extinction, 
making the intrinsic luminosity brighter by about 2 
%. 

Extinction is deduced from the distance between the 
observed color of the object and the intrinsic color 
of CTTSs on the color-color diagram. The intrinsic 
color of CTTSs deduced by [tje]cite.Meyer97Meyer et 
al.(1997) ([cite]cite.Meyer97Meyer et al.(1997)) has 
a dispersion as much as 0.17 mag in J — H, cor- 
responding to an uncertainty of 1.5 mag in visual 
extinction or 0.16 dex in bolometric luminosity. 
The intrinsic color of YSOs might be different from 
the intrinsic color of CTTSs. Central objects of the 
YSOs are M type stars whereas those of CTTSs arc 
mainly K type stars. Color of an M type dwarf is 
about 0.1 mag redder than that of a K type dwarf 
both in J - H and H - K. 

We assumed that the distance to the Taurus molecu- 
lar cloud is 140 pc. Hipparcos observations estimate 
142 pc with an uncertainty of ± 14 pc as the distance 
to the cloud (QDO [cite]cite.WichmannWichmann et 
al.(1998)). Therefore the uncertainty due to dis- 
tance is ± 0.1 dex in bolometric luminosity. 
The fractions of the photospheric luminosity are 
assumed to be 1.0, 0.8, and 0.5 of the total lu- 
minosity in the J-band. In these three cases, 
uncertainty is 0.3 dex in bolometric luminosity. 
If the fraction is smaller, photospheric luminos- 
ity is fainter. Although the fraction value could 
be deduced from veiling effect index rj (|(]00 
[cite]cite.Greene95Greene & Meyer(1995)), rj has a 
large uncertainty for all YSOs due to rather large 
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observational uncertainties in the equivalent widths 
of the metallic features. Therefore, we did not em- 
ploy rj to estimate the photospheric luminosity. 

4-. 4-2. HR diagrams 

Several groups have recently developed and refined 
evolutionary tracks for low-mass stars, brown dwarfs, and 
even giant planets. Figures [l(], [ll], and 12 show the HR 
diagrams of the YSOs and their companions. Overlaid 
are the evolutionary tracks of 0e]cite.DM94D'Antona 
& Mazzitelli(1994) ([cite]cite.DM94D Antona 

& Mazzitelli(1994)), §e]cite.DM98D'Antona & 
Mazzitelli(1998) ([citc]citc.DM98D'Antona & 

Mazzitelli(1998)), and |t]e]cite.BaraffeBaraffe et 
al.(1998) ([cite]cite.BaraffeBaraffe et al.(1998)). 

Mass and age for individual YSOs are estimated from 
the HR diagrams with the evolutionary tracks (table |^). 
Some YSOs are indeed very low-mass YSOs (typically 0.1 
Mq- 0.3 Mq). These objects are definitively a different 
population from CTTSs in term of their mass. The un- 
certainties in the mass are about a factor of 2 and those in 
the age are about a factor of 5 for the YSOs as shown in 
the figures. On the other hand, none of the mass and age 
of the companions are determined mainly due to the poor 
signal-to-noise ratio on the spectra. Additional notes are 
described in the appendix for most companions. 

Mass-luminosity relation of the YSOs is shown in 
figure [L3|, where also shown are the mass-luminosity 
relations of CTTSs (§00 [cite]cite.Strom88Strom 
et al.(1988)) and the predicted relations for 
three ages calculated from ^e]cite.BaraffcBaraffc et 
al.(1998) ([cite]cite.BaraffeBaraffc et al.(1998)). Most of 
the YSOs are distributed within one order of magnitude in 
the ages. Both CTTSs and the YSOs have near-infrared 
excesses, implying existence of circumstcllar disks. The 
ages of these objects indicate the survival time of the 
circumstellar disk, and it is interesting to ask whether 
the mass of the central object affects the evolution of 
circumstellar disk. The ages of most YSOs are less than 
10 7 yr, similar to the age of CTTSs. However, since 
the mass-age relation of YSOs depends strongly on 
evolutionary tracks (figures 0, [l|, and table |), we 
cannot conclude that disk survival time is independent 
on the mass of the central object. 

The small amount of objects prevents us from 
making a mass function of the YSOs. In the 
four other regions but except for Hcilcs Cloud 2, 
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Fig. 10. The HR diagram of YSOs with the 
evolutionary tracks of mc]citc.DM94D'Antona 

& Mazzitelli(1994) ([cite]cite.DM94D'Antona & 
Mazzitclli(1994)) with the Alexander opacities and the 
CM convection overlaid. Filled circles and filled squares show 
the loci of the YSOs deduced from the UKIRT observations 
and frorn the Subaru observations, respectively. T c ff(NIR) 
in table H are adopted. Error bars are tilted, because the 
bolometric luminosity is a function of th e J-band luminosity 
and the effective temperature (see §1.4.1. 



reach a mass limit of 0.03 A-/ Q for Av ~ with the survey. 
In order to take a census of brown dwarfs in Hcilcs cloud 
2 and to construct an initial mass function of the cloud, 
a deeper spectroscopic search is required. 

4-5. Comparison with other studies 

t'e]cite.WilkingWilking et 
al.(1999) ([cite]cite.WilkingWilking et al.(1999)) carried 
out if-band spectroscopic observations of low-luminosity 
sources in the p Oph cloud. They estimated spectral 
types of the objects using the Q index, then derived 
effective temperatures using the conversion for dwarfs. 
The effective temperature derived from their method is 
the same value as the lower limit of effective temperatures 
derived from the Q index in this paper. This is because 



§e]cite.Luhman00Luhman(2000) ([cite] cite. LuhmanOOLuhmaj^ 2§iQS&k have a higher temperature than M dwarfs even 
shows that the mass function has a peak with the same Q index, due to gravity dependence of the 
at 0.8 Af and is relatively flat between 0.1 water band depth (sec figure |). For example, ITG 40 has 
M Q and 0.8 M . |e]cite.Briceno98Briceno a Q of 0.67, which corresponds to a spectral type of M4.7 

(1998) ([cite]cite.Briceno98Briceno (1998)), as well as for Wilking's relation. This spectral type corresponds to 
|e]cite.Luhman00Luhman(2000) ([cite]cite.Luhman00Luhma&$gQ0$)in the dwarf scale, which is consistent with the 
suggest a deficit of brown dwarfs in the re- lower limit on the effective temperature derived from the 
although |^]cite. MartinOlMartin et 



gions, 

al.(2001) ([cite]cite.Martin01Martm et al.(2001)) have 
identified four brown dwarfs in two sparsely populated 
stellar groups. Given a limiting magnitude of the ITG 
survey (K = 13.4) and an age of 1 Myr, we expect to 



Q index in this paper. 

They estimated extinction of each object from the J — 
H color, assuming no infrared excesses in the J- and H- 
bands. This procedure could lead to a slightly higher value 
for extinction than our method. For example Av of ITG 40 



No. 



Taurus YSOs 

Table 4. Mass and age of the YSOs. 
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Ubject 


DM94' 


DM98 t 


BCAH98* 




mass 


age 


mass 


age 


mass 


age 




[Af ] 


[10 6 yr] 


[M ] 


[10 6 yr] 


[M Q ] 


[10 6 yr] 


ITG 2 


<0.15 


<2 


0.05 - 0.18 


<3 


<0.2 


<3 


ITG 6 


0.25 - 0.4 


20 - 50 § 


0.35 - 0.45 


>30 § 


0.4 - 0.5 


>50 § 


11 G 15A&U 


0.1 - 0.4 


0.1 - 3 


0.12 - 0.4 


0.2 - 5 


<0.55 


<8 


11G 15A 


0.1 - 0.3 


0.1 - 3 


0.12 - 0.3 


0.2 - 4 


<0.45 


<5 


lib 1/ 


<0.2 


<10 


0.08 - 0.3 


0.3 - 20 


<0.3 


<20 


lib ZT 


<0.15 


<2 


0.11 - 0.19 


1 - 5 


<0.2 


<4 


lib zt>A&;r> 


0.35 - 0.8 


1-10 


0.35 - 0.75 


2-8 


0.57 - 1.00 


4-18 


lib ZoA 


<0.35 


<1 


0.09 - 0.35 


<2 


<0.57 


<4 


ITG 25B 






>0.08 


>10 


>0.06 


>4.5 


ITG 27 


0.3 - 1.2 


1 - 5 


>0.3 


1 - 8 


0.57 - 1.2 


3-18 


ITG 33A&B 


>0.15 


>10 § 


>0.2 


>30 § 


>0.2 


>30 § 


ITG 33A 


>0.15 


>10 § 


>0.2 


>30 § 


>0.14 


>18 § 


ITG 40 


<0.2 


1 - 7 


0.14 - 0.3 


2-15 


0.1 - 0.32 


4-18 


ITG 41 














ITG 46 


>0.25 


>20 § 


>0.35 


>30 § 


>0.35 


>30 § 



[|=]cite. DM94D'Antona & Mazzitelli(1994) ([citc]cite.DM94D'Antona & Mazzitelli(1994)) with the Alexander opacitii 
invection 

b]cite.DM98D'Antona & Mazzitelli(1998) ([citc]cite.DM98D Antona & Mazzitelli(1998)) 
I 2]cite.BaraffeBaraffe et al.(1998) ([cite]cite.BaraffeBaraffe et al.(1998)) 
These ages are discussed further in the appendix. 
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Fig. 11. The HR diagram of YSOs with evo- 
lutionary tracks of ^e]cite.DM98D'Antona & 
Mazzitelli(1998) ([cite]cite.DM98D Antona & 
Mazzitelli(1998)) overlaid. 
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Fig. 12. The HR diagram of YSOs with evolutionary tracks 
of mcJcite.BaraffeBaraffe et al.(1998) ([cite]citc.BaraffeBaraffc 
et al.(1998)) overlaid. 
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Fig. 13. Mass - absolute J-band magnitude relation of the 
YSOs and the CTTSs listed in |e]cite.Strom88Strom et 
al.(1988) ([cite]cite.Strom88Strom et al.(1988)). The J-band 
magnitudes were corrected for interstellar extinction on the 
color-color diagram (ITG). The mass of the YSOs and 
the CTTSs are estimated from the evolutionary track of 
|t|e]citc.BaraffcBaraffe et al.(1998) ([cite]cite.BaraffeBaraffe et 
al.(1998)). The distance modulus of the Taurus molecular 
cloud is 5.73 (f]oO [cite]cite.EliasElias(1978)). 

is estimated to be 20.0 mag by their method. On the other 
hand, our value is 16.7 mag. Larger Av leads to higher 
values in bolomctric luminosity. Bolomctric luminosity 
of ITG 40 is estimated to be 0.1 Lq by their method, 
whereas 0.05 Lq by us, even though discrepancy is within 
uncertainty. 

tcJcitc.LuhmanRiekeLuhman & 
Rickc(1998) ([cite] cite. LuhmanRickcLuhman & 
Rickc(1998)) carried out i-T-band spectroscopy of 
low-luminosity sources in the L 1495 cloud in Taurus. 
They derived the spectral type of the object by fitting the 
2.21 fim feature and 2.26 /im strengths of the object to 
those of dwarf stars varying in amounts of veiling. With 
their method, spectral type of ITG 40 is estimated to be 
M3 or M7. Derived spectral type corresponds to 2600 K 
~ 3400 K, in agreement with our result. 

5. Conclusions 

1. We have carried out near- infrared spectroscopic ob- 
servations of 23 very low-luminosity YSO candidates 
and 5 of their companions in Heiles Cloud 2 in the 
Taurus molecular cloud. Near-infrared spectroscopy 
is essential to characterize objects in the color-color 
diagram. Out of the 28 objects, 5 objects have Br7 
in emission and 7 have "flat" spectra over the Br7 
feature. We conclude that these 12 objects are in- 



deed YSOs. 

2. Compiling near-infrared spectra of dwarfs and gi- 
ants taken by us as well as those in the literature, 
the ratio of the 2.26 fim feature to the 2.21 /im fea- 
ture turns out to be a good indicator of the effective 
temperature for M type stars. 

3. The effective temperatures of the YSOs are deter- 
mined. These objects are cool (T c g < 4000 K) 
YSOs. 

4. The mass of these YSOs is estimated from the HR 
diagram with recent evolutionary tracks. Some ob- 
jects are very low-mass (0.1 Mq- 0.3 Mq) YSOs. 

5. The age of these YSOs appears to be 10 5 - 10 7 years. 
However, the deduced age depends on evolutionary 
track models. 
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Appendix. Individual Objects 

ITG 2 

Mass of this object is estimated to be less than 0.2 
Mq by any evolutionary track, possibly a young brown 
dwarf. Age of this object is estimated to be less than 
3 x 10 6 . 

ITG 6 (GM Tau) 

' te] cite . GomezGomez et al.(1994) ([cite]cite.GomezGomez 
et al.(1994)) derived the magnitude of this star as 8.06 
in the .fT-band, one order of magnitude brighter than the 
magnitude derived by ITG and by |t]g]cite. KenyonKenyon 
& Hartmann(1995) ([cite] cite. KenyonKenyon & 
Hartmann(1995)). If the bolomctric luminosity is 
one order of magnitude brighter, the age of this star is 
estimated to be around 10 6 yr, while the mass of the star 
is st ill about 0.3 Mq. 

' te] cite . KenyonKenyon & 
Hartmann(1995) ([cite] cite. KenyonKenyon & 

Hartmann(1995)) classified spectral type of this star 
as "continuum". Therefore the spectrum of this object 
has a large amount of veiling. Since only the lower limit 
of the effective temperature is estimated from the Q 
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index, true effective temperature might be lower than 
that estimated here. If half of the if -band flux comes 
from a circumstcllar disk, for example, the Q index 
changes from 0.81 to 0.62 and the effective temperature 
to 2700 K ~ 3500 K. This value is not inconsistent 
with the effective temperature deduced from the ratio of 
the 2.26 /im feature to the 2.21 /im feature. With this 
temperature, the mass and the age of the object would 
be 0.1 M n and 3 x 10 6 



yr. 



ITG 9 



The fact that both spectra of the binary have neither 
the strong metallic features nor the deep water absorption 
bands indicates high effective temperatures of the objects. 
They both might be field-stars. 

ITG 15 A 

The effective temperatures deduced from the UKIRT 
observations and that from the Subaru observations are 
very consistent with each other. 

ITG 25 (IRAS 04370+2559) 

While the UKIRT observations and the Subaru obser- 
vations are consistent with each other in the strength of 
the 2.26 /mi feature as well as the Q index, there is a 
discrepancy in the strength of the 2.21 /im feature, result- 
ing in a different effective temperature. Because of higher 
signal-to-noise ratio in the Subaru observations than in 
the UKIRT observations, we use the effective temperature 
derived from the Subaru observations. For the Subaru ob- 
servations, effective temperatures derived from the feature 
ratio and from the Q index are also different. As is the 
case of ITG 6, the veiling effect can easily increase the Q 
index, leading to a higher effective temperature. Because 
this object is an IRAS source indicating its young age, we 
believe the effective temperature derived from the feature 
ratio rather than that from the Q index. 

For consistency between the effective temperature de- 
rived from the feature ratio and that from the Q index, 
the Q index should be 0.65, which means that 70% of the 
isT-band luminosity is not of photosphcric origin but from 
a circumstcllar disk. If the same portion of the luminosity 
is due to the veiling for the companion ITG 25B, the Q 
index of the companion is suppressed to 0.08. This very 
low value indicates that effective temperature is less than 
2700 K and the age is less than 1 x 10 6 yr, implying a 
low-mass brown dwarf. 

ITG 33 



of HK Tau (§00 [cite]cite.Koresko98Koresko(1998); §00 
[cite]cite.Stapelfeldt98Stapelfeldt et al.(1998)). If so, 
visual extinction would be very large as much as 50 mag, 
like the HK Tau companion, inconsistent with A v for 
ITG 33A. The location of ITG 33A is odd as a YSO in 
the HR diagram. Veiling effect may lead to high effective 
temperature. A spectrum with high signal-to-noise 
ratio is required to estimate effective temperature from 
met allic features. 

te]cite.Martin00Martm(2000) ([cite]cite.Martin00Martm(2000)) 
also classified ITG 33B as a field-star from optical spec- 
troscopy. However, their optical spectrum, as well as 
our near-infrared spectrum, are too noisy to classify this 
object as a field-star or a featureless YSO. Higher signal- 
to-noise ratio spectrum at optical and/or of near-infrared 
wavelength is required. 

ITG 41 (IRAS 04385+2550) 

This object is identified at optical wavelengths as Haro 
6-34. Weak absorption features which may imply large 
amount of veiling prevents us from estimating the effective 
temperature, thus its mass and age. 

ITG 45 

c ^e]cite.Martin00Martm(2000) ([citc]cite.Martin00Martm(2000)) 
found an Ha absorption line in a spectrum of ITG 45A, 
indicating a field-star. Near-infrared spectrum of the pri- 
mary contains Br7 absorption, also indicating a field-star. 
@e]cite.Martin00Martfn(2000) ([cite]cite.Martin00Martm(2000)) 
claims that the companion, ITG 45B, is also a field-star. 
However, the optical spectrum, as well as our near- 
infrared spectrum, is too noisy to classify it as a field-star 
or a YSO. 

ITG 46 

The location of ITG 46 in the HR diagram is unusual. 
Again a large mount of veiling may lead to high effective 
temperature. Otherwise, this object may be a field-star. 
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